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Soliton interaction in periodically twisted birefringent optical fibers has 

been analysed analytically with refernce to soliton switching. For this purpose 

we construct the exact general two-soliton solution of the associated coupled 

system and investigate its asymptotic behaviour. Using the results of our 

analytical approach we point out that the interaction can be used as a switch 

to suppress or to enhance soliton switching dynamics, if one injects multi- 

soliton as an input pulse in the periodically twisted birefringent fiber. 

PACS Numbers 42.81Dp, 42.65Tg, 03.40kf 

It is a well-known fact that single-core fiber supports two distinct modes of propagation 

as a result of birefringence effect, which can be introduced through twisting the fiber during 

the preform stage of formation or through stress induced birefringence mechanism. For 

many years the propagation of solitons in stress induced birefringent nonlinear Kerr media 

with reference to optical fibers is an intensive research subject for theoretical as well as 

experimental investigations [1]. The topic of propagation in twisted birefringent optical 

fiber is also gaining considerable interest theoretically and experimentally in recent times as 

most of the nonlinear directional couplers are based on such fibers. One may mention the 

example of the rocking rotator [1] which can be used as a switching device at high-power [2] 

and as a filter at low-power [3]. 

Soliton propagation in the periodically twisted birefringent fiber is usually described by 

using the coupled nonlinear Schrodinger (CNLS) family of equations [1,2]. However such 

CNLS equations are in general not completely integrable. Interestingly if we assume the 
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value of ellipticity angle to be 35° then the dynamics of soliton interaction in a periodically 
twisted birefringent fiber coupler can be described by coupled wave equations of the form 

[4] 



iqiz + qitt + pqi + m2 + 2Ai( gi ^ + 


1 ?2 n^i 


= 0, 




n2z + ?2« - p?2 + mi + 2/i( 1 1^ + 


1 ?2 ?)q2 


= 0, 


(1) 



where qi (z,t) and 52 are slowly varying envelopes of two orthogonally polarized modes, 
z and t are respectively the normalized distance and time and x and p are the normalized 
linear coupling constants caused by the periodic twist of the birefringence axes and the 
phase-velocity mismatch from resonance respectively. If the linear coupling constants are 
absent (that is, x — p — Q), then one can easily recognize the system(l) to be the celebrated 
integrable Manakov model [5]. The resulting Manakov equation is receving renewed atten- 
tion recently as it describes the effects of averaged random birefringence on an orthogonally 
polarized pulse in a real fiber [6] . When the birefringene axes of fiber are periodically twisted 
during the drawing process there is a periodic intensity exchange between the orthogonally 
polarized modes [1-4] and it can be modelled by Eq.(l) if the value of eUipticity angle is 35° 
[4]. The linear coupling length where the maximum power is transfered from one mode to 
the other is tt 

/(2^x2 + p2). At the resonance wavelength, the linear parameter p = and 
the linear coupling length increases to tt/ (2x). The schematic of the experimental apparatus 
used to observe switching characters in the periodically twisted birefringent fiber is given in 
Ref. [2]. Particularly the dependence of switching characters on the input power, operating 
wavelength, twist magnitude and twist period are described for example in Refs. [2] and [4] . 

Considerable attention has been paid in the literature [1,5,7,8] to study soliton collision in 
the birefringent fiber. Particularly by using the system (1) in the absence of linear coupling 
terms Manakov [5] pointed out that during soliton collision their velocities and amplitudes 
(intensities) do not change but the associated unit polarization vectors do change provided 
they are neither parallel nor orthogonal (see also the paper of Silmon-Clyde and Elgin [8] 
for a discussion in terms of Stokes vectors). Further Menyuk [7] has observed at the value 
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of the ellipticity angle 35°, where the Manakov equation holds good, that a soliton of one 
polarization when interacting with a switching pulse of the other polarization does not 
develop a shadow and also does not change shape. However very recently we have proved [9] 
by constructing the most general two-sohton solution of the Manakov model that it has the 
property that soliton in birefringent fiber can in general change its shape after interaction 
due to change in intensity distribution among the modes enenthough the total energy is 
conserved. In this letter we investigate the implication of this property of the solitons when 
the additional effects due to the periodic rotation of birefringence axes are included by 
constructing the exact two-soliton solution of the system (1). In particular we point out 
that interaction can be used as a switch to suppress or to induce soliton switching, if we 
inject multisoliton as an input pulse in the periodically twisted birefringent fiber. 
The coupled system (1) reduces to the celebrated integrable Manakov model [5], 



iqiMz + qiMtt + 2At( qiM + 


q2M n?iM 


= 0, 




iq2Mz + q2Mtt + 2At( 1 gfiM r + 


1 |2\ 

1 ?2M 1 )q2M 


= 0, 


(2) 



(the subscript M refers to the Manakov model) under the transformation [4] , 

qi = cos(^/2)e'r^giM - sin(^/2)e-^^^g2M, 

q2 = sin(^/2)e*r^giM + cos(^/2)e-*r^g2M, (3) 

where T — (p^ + n^yl'^ and Q — tan~'^(x/p). Belanger and Pare [10] and also briefly in [1] 
have shown that the system (1) without linear self-coupling (p = 0) has simple solitary wave 
solutions exhibiting energy exchange between the modes. By using the one-soliton solution 
of the Manakov model (2) in (3), Potasek [4] has pointed out the possibility of periodic 
intensity exchange between the orthogonally polarised modes qx and q-2, in the coupled system 
(1), when both the linear coupling constants x and p are present. An interesting question 
arises here when one considers multi-soliton solutions for the Manakov model that admits 
both elastic and inelastic (shape changing) type of collisions depending upon the initial 
conditions or arbitrary parameters as shown in Ref. [9]. Then how does the switching and 
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energy sharing properties get modified for the multi-sohton solutions of the system (1)? We 
show in this letter that indeed novel features in the intensity sharing and different switching 
properties do arise when the most general two-soliton solution is considered for (1). 

For our analysis we make use of the general two-sohton solution of the Manakov system 
(2), reported in [9], and obtain the corresponding two-soliton solution of (1) through (3). It 
has the form : 

qi = {[cos{e/2)e'^'ai - sm{e/2)e-'^'p,]e'^' + [cos{e/2)e^^'a2 - sin(^/2) 
e-''^'p2]e'^^ + [cos(^/2)e^^^+''i - sm{e/2)e-''^'+^'^]e'^'^''*^+'''+ 
[cos(^/2)e'^^+'^2 - sin(^/2)e-^^^+''2]e'?i+'?2+'?2}/D„, 

q2 — as qi above with the replacements cos(^/2) sm{9/2) and 

sin(^/2) ^ -cos(^/2), (4) 

(the symbol * denotes complex conjugation). Here Dn — I + exp(77i + rj^ + Ri) + exp{r]i + 
V2 + ^o) + exp(r/]' + r/2 + 5^) + exp(r/2+ rj^ + R2) + exp(r/i + + V2 + V2 + -^s) and rjj = 
kj(t + ikjz), j = l,2. The parameters exp(5i) = (fci — k2){aiK.2i — Q;2fi;ii)/(A;i + k^) {k^ + 
k2),exp{S2) = {k2 - ki){a2Ki2 - q;iK22)/(^2 + k^){ki + A;^),exp(5^) = {ki - /C2) {PiK'2i - 
(^2i^ii)/{ki + kl){kl + A;2),exp(52) = (^2 - A;i)(/32«i2 - (^ii^22)/{k2 + kl){k^+ k^), exp(5o) = 
Ku/{ki + k2),exp{Ri) = Kn/{ki + kl),exp{R2) = K22/{k2+ fca); exp(i?3) =| ki - k2 ^ 
(«ii«22 - «i2«2i)/(^i + kl){k2 + kl) \ki + kl p and = //(ctjO;* + l3i/3*){ki + i,j = 

1, 2. The six arbitrary complex parameters cci^ 0:2, (32, fciand ^2 determine the amplitude, 
velocity and phase of the asymptotic soliton forms of (4). Now in order to bring out the 
nature of the solitons of (1) and their interactions including exchange of energy we carryout 
an asymptotic analysis of the solution (4). To be specific we choose the arbitrary complex 
parameters ki,i = 1,2, as ku > k2j, kiR > and k2R > (here subscripts I and R refer to 
the imaginary and real parts). 

I. Limit z— > —00 : As —00 we can identify two independent solitons denoted by 
soliton 1 and soliton 2 with the above choices of A;iand ^2- Soliton 1 will be centered 
around tjir — kiR{t — 2kuz) ~ (when r]2R —00) and soliton 2 will be centered around 
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r]2R = hnit - 2k2iz) ~ (when r]iR oo) . 

(a) Soliton 1 {tjir ~ 0, r]2R -oo) : 

q, - [cos{e/2)e'^^Alj,-sm(e/2)e-'^^Al-^]q'-, 

= [sm{e/2)e^^^ Al-^ + cos{e/2)e-'^^Al-^]q'- , (5) 

where q^' = kiRexp{ir]ii) sech{r]iR + Ri/2),r]u = kut + {kfj^ - kfj)z, (A}^^ Al]^) = [fj,{ 
q;iQ!i + /?i/3i)]"^/^(q;i,/3i) and | A\]^ ^ + | A^]^ p= If/j, . Here ^2m) refers polarization 

unit vector of Manakov one-sohton solution, which can be obtained from (5) by substituting 
^ = r = 0, superscripts 1- denote soliton 1 at the limit — oo and subscripts 1 and 2 refer 
to the modes giandg2- Eq.(5) exhibits the same form of one-soliton solution of (1) reported 
by Potasek in the Ref. [4]. If we parametrize [8] the unit polarization vector as {Alj^^ 
^2m) — (*^o^(^p ) cxp(iap]^), sin(^^p^) cxp(iap2 )) then we can identify 9p~ as the polarization 
angle and the phases a^J^ ^ ot^^ corresponds to the state of elliptical polarization. 

(b) Soliton 2 {r]2R ~ 0, r]xR oo) : 

q, - [cos(^/2)e^r^A2- -sin(^/2)e-r^A2-]g2-, 

q^ = [sin(^/2)e*r^A2- + cos{e/2)e-'^^Al-^]q'- , (6) 

where = k2Rex.p{ir]2i) sech{ri2R + {R3 - Ri)/2), r]2i = k2it + - k^j)z, {A^^^ Al'^) = 
(ai/a*) c[ii{a2a* + /32f3;)]-^/^[{ai^(3i)K^^ - {02,^2)^21], \ Allj p + | Al]^ p = in which 
ai = {ki + k2)[{ki_k2){ala2 + ^1^2)]^^^ and c=[ 1/| K12 ^ - l//«ii/«22]^^^- It is interesting 
to note from (5) and (6) that the form of A^'^ and Al'^ in (6) differs from the values Al]^ 
and Al]^ in (5) and the former contains more number of parameters, eventhough (6) is an 
exact one-soliton solution of the system (1) just like the solution (5). Further in the special 
case ai : a2 — Pi '■ P2, the form of (6) reduces to the form of (5) with parameters specifying 
soliton 2. So Eq.(6) may be considered as the most general one-soliton solution form of (1). 

II. Limit z— > 00 : We now analyse the form of the solitons after interactions as z— > 00 : 
(a) Soliton 1 (771^ ~ 0, ri2R 00) : 

gi - [cos{e/2)e^^^A\+, -sm{e/2)e-^^^Alt,]q'+, 
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?2 = [^9/2)6^^^ Alt, + cos{9/2)e-^^^Alt,]q'^, (7) 

where = ki^exp(iriu) sec h(r]iR + (R^ - i?2)/2), (A}^ ) = (aa/a^) c[//(q;iq;^ + 

_ (a2,/32)K22^],| A}+ p + | p = in which 8.2 = {k^ + 

A;i)[(A;i_A;2)(q;iQ!2 + /3i/32)]^^^- Note that "^^im A2][^ 7^ ^\mi except when ai : 
0^2 — Pi '■ P2, corresponding to pure elastic coUision in the Manakov model [9]. 
(b) Soliton 2 {r]2R ~ 0,riiR -00) : 

gi = [cos(^/2)e^r^A2+ -sin(^/2)e-rM2+]g2+, 

52 = [sm{e/2)e^^^ All, + cos{e/2)e-^^^All,]q'+, (8) 

where (f+ = /c2J?exp(i772/) sec h{rj2R + R2/2), {Afl^ AH^) = [//(a2a;2 + /92/?2)]"^^^(«2,/92) and 
I ^iM P + I ^2M P= Here also 7^A^~ and 7^ unless ai : q;2 = /9i : /32- 

Now we recognise from (5-8) that not only the phase-factors but also the overall shapes 
of the solitons get modified (due to the intensity redistribution among the solitons) after 
undergoing interaction when the two coupled one-solitons (5) and (6) move from z— > —00 to 
z— s> cxo as shown in Eq.(7) and Eq.(8). To facilitate the understanding of the above behaviour 
with reference to the optical soliton switching between the orthogonally polarised modes, it 
is convenient to obtain the oscillating parts of the intensities associated with the asymptotic 
forms (5-8) as 

' = I A[S r cos2(^/2)+ I r sin^(^/2) + (-1)' I AS II ^Iti I 

sin(^) cos(2rz + 0"^), l,j=l,2 (1 7^ j), z ^ ipoo, (9) 

where 0"^ = tan~^ (^^) — tan~^ (^f^) ■ The presence or absence of the last term involving 
the factor cos(2r2; -|- (fP^) plays a crucial role in the switching behaviour of solitons as 
demonstrated below. 

As we have mentioned before the values of A^j^ (j,n=l,2) change to new values A^j^ 
due to the collision between two co-propagating solitons, namely soliton 1 and soliton 
2, without violating the condition | p -|- | A^^j P= The amount of change 
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q^^{z,t) 



(A"jj^— A"^) can be estimated by assigning suitable values to the arbitrary parameters 
ki,k2,Q;i,Q;2,/3i, and /32, appearing in the expressions for (j,n=l,2). Further from (9) 
it is obvious that depending upon the values of |A^^ |, the nature of switching dynamics 
supported by the system (1) also changes. The change in can be parmeterized as 
(cos(6'p^) exp(iQ;pj^), sin(6'p^) exp(iQ;p^)), where as long as the phases a^^ ^ a^^ the state 
of polarization is preserved during interaction. Therefore in general without affecting the 
state of polarization, the switching dynamics can be changed just by changing 74"^with the 
help of the polarization angle. In the following we briefly discuss the different changes which 
can occur in the intensity exchange between qi and q2 modes with respect to soliton 1 and 
soliton 2 due to the above mentioned coUision by considering |^| and |^| defined in (9) 
for each 1=1,2 and n=l,2 values. 

Case 1 : All the | |'s (j,n=l,2) are nonzero. In this case due to the presence of 

cos(2r2; + 0"'^) term on the right hand side of (9), there is a periodic intensity switching 
which is always present in both the solitons and in both the components before as well as 
after the interaction. Of course the conservation relations | A^^ P + I P~I P 

+ I >l2M for the total intensity are always vahd. However the switching dynamics 

appearing before and after interaction is not similar in form, due to the condition A'^'^ 
y^Ajf^, j,n=l,2 , except when ai : a2 = Pi P2 as mentioned before, giving rise to a partial 
suppression or enhancement of the periodically varying intensities. 

Case 2 : Any one of the | |'s is zero and others are nonzero. For example if | A\'^ |~ 
0, corresponding to the condition 0:2(^2+^2) ('^i'^2 + = C(;2 P + | P2 ^(^i + ^2)) 

then the switching in the intensity of the soliton 1 gets fully suppressed in both the modes 
qi and 52, while it persists for the other soliton. This is illustrated in Fig. la for the choosen 
parameters values, namely ki = 1 + i, k2 = 2 — i, ai = Pi = (32 = ^, (^2 = (39 + i80)/89(~ 
exp(i64°)),p = 0.25 and x = 0.5, for which | Al'^ |~ 0.7,| |~ 0.5, | A^'^ |~ 0.7,| Al]^ |~ 
0.86,1 Al-j^ |~ 0.06,1 Afl^ |~ 0.7, | A^-^ |~ 0.99 and | ^2M I~ 0-7, satisfying the condtion 
that I AI'Ij |~ 0. Similar phenomenon can be seen if | A^'lj |~ instead of | |~ 0. But 
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if we choose | A^~l^ |~ (j=l or j—2) the periodic intensity exchange with respect to sohton 
2 will be suppressed, while it persists in soliton 1. On the other hand, if | Aj^ |~ (j=l or 
i—2; n=l or n=2), then there is no switching in the intensity of soliton n before interaction, 
but the switching appears after interaction in that soliton and so there is an inducement of 
switching due to the inteaction. Thus the interaction itself acts as a switch to suppress or 
to enhance the switching dynamics. 

Case 3 : Any two of the | ^4"^ |'s are zero and others are nonzero without violating 
the conservation conditions. For concreteness, let | A^^ \ ~ (n = 1,2), which imphes 
the condition cti ~ q;2 ~ 0. It implies that | \ (n=l,2) should also simultaneously 

vanish. Consequently there is no switching between the modes qi and q2 either before or 
after interaction and there will be only inelastic (shape changing) scattering as discussed for 
the Manakov model in Ref. [9]. Similar observations can also be made if one makes any two 
of the I Aj^ I 's to vanish. However, one can identify the interesting possibility of switching 
existing in soliton 1 only before interaction, which gets interchanged with soliton 2 after 
interaction, by allowing one of the two | \ 's and another one of | A^'l^ \ 's simultaneously 
to take the value zero, corresponding to the condition [| ki + k2 P (| ai p + | Pi p)(| a2 P 
+ I P2 ni = [(^1 + k\){k2 + k2) I aio.^ + /3i/92 P]- Figure lb, for the choosen parameters 
namely /ci = 1 + iO.l, /cs = 1 - iO.l, ai = 0.86 + iO.5, ^2 = 0.5 + iO.86, pi = 0.7 + iO.72, 
P2 = 0.44 + iO.9, and p = X = 0.25, for which | A\]^ |~ 0.7, | Al]^ |~ 0.05, | Al]^ |~ 0.7, | 
^2M h 0.99,1 Al+^ |~ 0.04,1 |~ 0.7, | Al+^ |~ 0.99 and | Al+^ |~ 0.7 satisfying the 
above condtion, shows that the interaction induces the periodic intensity exchange between 
the two modes of soliton 2 while it suppresses the switching dynamics in soliton 1. 

To conclude, by studying the interaction between two co-propagating solitons in the 
periodically twisted birefringent fiber with reference to soliton switching, we have observed 
several possible ways to use interaction as a switch to suppress or to induce the switching 
dynamics. The basic underlying mechanism of such possibilities is the inelastic (shape 
changing) nature of the soliton interaction which arises essentially due to changes in the 



8 



polarization angle and so in the overall amplitude of the solitons. Since MoUenauer et al [8] 
have demonstrated polarization scattering by soliton-soliton collision, it shoud also possible 
to experimentally study the phenomena described in this paper by using specially fabricated 
optical fibers. These possibilities should have important ramifications in nonlinear switching 
devices like the rocking rotator. 
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Part I 

Figure Caption 

Fig.l : Typical evolution of the intensity profiles | gi P and | q2 P of the two-soliton solution 
(4) (a) showing the suppression in switching between the two modes of si soliton and (b) 
showing the suppression in switching of si soliton and enhancement in s2 soliton (while 
undergoing a large phase shift) for the parameter values given in the text. 
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